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Abstract Atomic oxygen density in the upper thermosphere (~300 km) can be calculated 
using around based incoherent scatter radar and Fabry-Perot interferometer measuremeiits. 
Burnside et al. [1991] was the first to try this method, but Buonsanto et al. provided an extensive 
treatment of the method in 1 997. This paper further examines the method using 46 nights of data 
collected over six years and the latest information on the oxygen collision fi^equency. The 
method is compared with the MSlS-86 atomic oxygen prediction values, which are based upon 
in situ rocket bom and satellite measurements from the 70's to the mid-80's. In general, the 
method supports the MSIS-86 model, but indicates several areas of discrepancy. Furthermore, 
no direct correlation is found between the geomagnetic conditions and the difference between the 
method and MSIS-86 predictions. 


Introduction 

Due to its extreme height, elusiveness has often pervaded study of the thermosphere. In 
this respect, measurement of neutral atomic oxygen density ([O]) has been no different. Because 
of the lack of a constant data supply, researchers have had to rely upon the MSIS-86 model tor 
calculations requiring thermospheric oxygen densities. While Hedin's model has proven to be 
generally reliable, it has limits and significant error margins. As indicated by Hedin [1987] and 
Grossmann et al. [2000], MSIS predictions below 200 km become increasingly divergent from 
realitv due to the lack of data for this region. Grossmann's results from the CRISTA experiments 
indicated an average of 40% less [O] than predicted during low solar activity. Buonsanto et al. 
[1997] estimated a 20% uncertainty for MSIS predictions higher in the thermosphere based on 
the information in Hedin [1987]. Burnside et al. [1991] found large errors during 
geomagnetically disturbed periods above Arecibo at 300 km with measured densities twice the 
predictions. 

Such uncenainties for [O] predictions are unsatisfactory because atomic oxygen plays an 
important role in various thermospheric reactions. At an altitude of 120 km, neutral atomic 
oxygen constitutes roughly 20% of the total density of the atmosphere. Here it affects the heat 
transfer between the mesosphere and low'er thermosphere by recombining into O 2 and 
deactivating CO; [Ward et al., 1993]. At an altitude of 300 km, neutral atomic oxygen 
constitutes 'roughly 80% of the total density of the atmosphere. Here it moves solar energy 
deeper into earth's theimosphere by transferring an electron through strong resonance charge 
exchange interaction with O', providing the principal method of heating the thermosphere 
[Omidy ar et al., 1998]. 

Within the last decade, several methods have been introduced to calculate neutral atomic 
oxygen density for various areas of the thermosphere, including twilight 7320 A airglow 
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emissions [McDade et al., 1991], satellite far ultraviolet images (DE 2) [Drob et al., 1999], 
satellite 63 fim infrared measurements (CRISTA) [Grossmann et al., 2000], and a host of rocket 
bom measurements for the 80 to 200 km range [partial list: Shepherd et al., 1995; Kita et al., 
1996; Melo et al., 1996]. A new technique for determining [O] at about 300 km using ground 
based incoherent scatter radar (ISR) and Fabry-Perot interferometer (FPI) measurements was 
introduced at the beginning of the decade. Burnside et al. [1991] was the first to use this 
technique, which he applied to data collected at Arecibo. Buonsanto et al. [1992] provided a 
fuller treatment of the technique at Millstone Hill for fourteen nights of data ranging from March 
19, 1990 to May 16, 1991. In this paper, we will further examine this technique for 46 nights 
with a larger data set, a greater span of time, and the most recent information on O -O collision 
frequency. 

Data Collection and Instrumentation 

The data selected provides a wide spectrum of solar and geomagnetic activity. The data 
set includes 46 nights and 899 points of coincident ISR and FPI measurements from Aug. 19, 
1988 to Mar. 3, 1995 at Millstone Hill Observatory in Massachusetts, USA (42.6° N and 71.5° 
W). Table 1 lists every night of collected data. Buonsanto et al. [1997] originally processed the 
data. 

By employing Thomson backscatter from ionospheric electrons, the incoherent scatter 
radar (ISR) is able to ascertain ion velocity along the magnetic field (V||), ion (electron) density 
(Ni), ion temperature (Ti), and electron temperature (T«). The Fabry-Perot interferometer (FPI) 
observes the atomic oxygen at 630 nm and calculates the neutral wind velocity (U) by reading 
the Doppler shift of the nightglow. Although the FPI can also provide neutral temperature 
measurements, these tend to have large uncertainties, making them less reliable than the radar 
measurements. The ISR is able to operate nearly continuously, but poor visibility from cloudy 
nights limits data from the FPI. 

By combining the ISR electron density profile with data from the MSIS-86 model, the 
altitude of the maximum 630 nm emissions was estimated for each data point [Link & Cogger, 
1988]. It is at this altitude that the FPI makes its neutral wind velocity (U) readings, so all ISR 
data was selected for this altitude. Furthermore, the MSIS-86 atomic oxygen density predictions 
were also based upon this altitude. The maximum emission altitude tends to range about 50 km 
on any one night and may be as low as 250 km or as high as 350 km. It is generally about 300 
km. 

For more information on the data collection and instrumentation, we recommend Sipler et 
al. [1991] and Buonsanto et al. [1997]. Sipler provides a complete discussion of the FPI and a 
summary of the ISR, and both Sipler and Buonsanto nicely overview the data collection for 
substantial portions of the same data used in this study. 

Theory: From Diffusion Velocity to Oxygen Density 

The diffusion velocity of the ions (Vd) can be calculated through two different forms; 

V^=u cos I — V[| (1) 
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where Tp = (Tj + Te) / 2 is the plasma temperature, U the horizontal meridional neutral wind 
velocity, I the magnetic dip angle, Vy the ion velocity along the magnetic field, k the Boltzmann 
constant, rrti the ion mass, Vi„ the ion-neutral collision frequency, Ni the ion (electron) density, Z 
the altitude, and g the acceleration due to gravity [Reddy et ah, 1994]. It is assumed that the 
ionosphere is composed of O^, so all calculations refer to this ion. 

The collision frequency is an important component of Equation (2). At 300 km, the main 
contributions to Vm are 0^-0, 0^-N2, and 0^-02 collisions, of which the 0^-0 collision 
frequency is the largest because of the high concentration of neutral atomic oxygen. Due to the 
difficult nature of experimentally determining Vm at low energies, we must use a theoretical 
formula instead; 
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< = Qf[0] + 6.9 X 10-'^[iV2] + 6.7 X 10-'^[O] , <4) 

where Qd is the theoretical collision cross section of 0^-0. Because N2 and O2 are minor 
constituents of the thermosphere at 300 km (roughly 19% aAd 1%), their collision frequencies 
are approximated using Banks [1996] and MSIS-86 density predictions without appreciable 
inaccuracies. 

Conversely, the V0+.0 has a large impact on the results. Its exact value has been the topic 
of much research and debate, especially over the last decade. For a considerable time, theoretical 
results varied significantly, but the recent correlation in the formulas by Stubbe [1968], Stallcop 
et al. [1991], and Pesnell et al. [1993], has provided a firmer foundation. The formula by Pesnell 
et al. is given by: 


=3.0xl0-'’r/"(l-0.1351og7;)' ( 5 ) 

where Tr = (Ti + T„) / 2, with Tn the neutral temperature and T3 = Tr / 1000. Pesnell's formula is 
larger than the commonly used formula of Dalgamo [1964] by 1.38, 1.32, 1.28, 1.24, and 1.22 at 
500 K, 750 K, 1000 K, 1250 K, and 1500 K respectively. 

The quest for empirical verification of Qa started roughly, but has began to come into 
focus. Burnside et al. [1987] used ISR and FPI data from Arecibo to calculate Vo+.q, but his 
values were 1.7 times larger than those predicted by Dalgamo. Sipler et al. [1991] used a similar 
technique at Millstone Hill to calculate values 1.9 times larger than Dalgamo, and Salah [1993] 
recommended an official correction value of F = 1.7, which was termed the Burnside Factor. 

In the next few years, a series of papers found fault with the previous data analysis 
methods. Reddy et al. [1994] demonstrated a systematic overestimate of Q<j from errors in the 
measurements. Hines et al. [1997] questioned the use of the least-squares method in this 
application. Buonsanto et al. [1997] conducted Monte Carlo simulations on Millstone Hill data, 
discovered a systematic overestimate from random errors, and calculated an unbiased F = 
1.4±0.3. Finally, Omidvar et al. [1998] conducted three different analysis methods on the 
Millstone Hill data and determined it supported the theoretical formula by Pesnell et al. 
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If the collision frequency is now assumed to be known, we can determine neutral atomic 
oxygen density in the thermosphere using Equations (1) and (2) as follows: 
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O is the neutral atomic oxygen density, J groups many of the terms from Equation (2), Vcorrection 
accounts for the Vq+-n 2 and Vch.o 2 , W is the diffusion velocity from both ISR and FPI wind data, 
and D is the theoretical portion from ISR temperature data. 

Analysis 

The past two studies of atomic oxygen density by Burnside et al. [1991] and Buonsanto et 
al. [1992] used a correction factor F = 1.7 and Dalgamo's formula for Qd. Based upon Omidvar 
et al. [1998], the present work uses a correction of F = 1.0 and Pesnell's formula for Qd. We 
have adjusted all of Buonsanto's data accordingly. 

To provide a method of evaluating the data with itself, we divide the original 899 points 
into two data sets: Set A from August 19, 1988 to March 23, 1990 with 445 points and Set B 
from March 29, 1990 to March 3, 1995 with 454 points. After removing all outliers that produce 
F values below 0.35 or above 3.0, Group A has 398 points and Group B has 390 points. Each 
point consists of a D, W, and O. The D and W are derived from empirical data, while the O is 
the MSIS-86 [O] prediction for the particular time, position, and altitude of the matching data. 
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In a perfect world, the MS IS derived [O] will equal the DAV derived [O], so we compare these 
two values to check the accuracy of using ISR-FPI data to measure [O], check the validity of the 
MSIS model, and check the accuracy of the theoretical value of Qd- All comparison is 
conducted at the 5% level of significance. 

The standard deviations for both data Sets A and B (Table 2) demonstrate their wide data 
spreads (see Figure 1 and 2). Both data sets are comparable. As can be seen in Figures 3 and 4, 
the DAV's from Set A and Set B are heavily skewed toward the right side, approximating a 
lognormal distribution. In order to perform correlation analysis upon the data, which requires a 
Gaussian distribution, we take the log of all the D/W data points. This transforms the data from 
an approximate lognormal distribution to an approximate Gaussian distribution (see Figures 5 
and 6). The MSIS-86 O predictions already form Gaussian distributions. 

The correlation analysis comparing O with DAV returned r^ = 0.237 and r^ = 0.139 for 
Set A and Set B respectively, indicating a low, but significant, association in their variances. In 
general, though, D/W does not follow the trends of O. The correlation coefficients produced 
from each night of data range from r = 0.93 to r = -0.80. Ten nights display a significant positive 
correlation. One night displayed a significant negative correlation. Twenty-four nights have no 
significant correlation. The remaining ten nights have less than six data points, so we did not 
calculate their individual correlation values. 

Although O and D/W are of similar magnitudes, many data points have substantial 
differences. In Set A, D/W is at least ±25% greater than O for 217 points. In Set B, D/W is at 
least ±25% greater than O for 148 points. 

Perhaps the best analysis of the data is on a night by night basis. Twenty-four nights 
have average percent differences of the absolute value between D/W and O greater than 40%. 
Nonetheless, data from eight nights closely follows MSIS predictions. One of the best examples 
is September 12, 1991, which has an average difference of 21.2% and a significant correlation of 
r = 0.503 (Figure 7). 

For the majority of the nights, the data roughly approximates the MSIS-86 predictions but 
includes large deviations. On a few of these nights, the data is extremely erratic, with December 
29, 1989 being one of the worst examples (Figure 8). The erratic nature of the density on this 
night results from both dramatic shifts in wind speed and large changes in D. A number of 
nights have consistently different trends from the MSIS-86 predictions, such as in Figures 9, 10, 
11, and 12. 

We found no significant positive correlation between the geomagnetic activity and the 
average percent difference between D/W and O (Figure 13). For example, April 10, 1990 
exhibits increased densities during a geomagnetic storm reaching a max Kp of 8- (Figure 14), 
but so does January 26, 1995 even though it maintains a low Kp of 1 (Figure 15). This is 
somewhat in contrast to Burnside [1991] who suggested densities dramatically increased over 
MSIS-86 predictions during geomagnetically disturbed periods based upon measurements on 
July 14, 1985 (max Kp 4+) and January 15, 1988 (max Kp 7+). 

Over half of the data was collected from Aug. 19, 1988 to September 21, 1990. This 
concentration of twenty-six nights over two years allows the opportunity to search for seasonal 
trends. We average each night of data to approximate the atomic oxygen density for that period 
of the year. Figure 16 shows a clear seasonal cycle with a winter maximum and summer 
minimum. In general, the alternative method produced averages larger than the MSIS-86 
predicted averages. 

Sources of Error 
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Error in the measurements of the two wind velocities (ion velocity from the ISR and the 
neutral wind velocity from the FPI) will have the largest effect upon the results. Both of these 
velocities factor into the calculation of the diffusion velocity (W). Buonsanto [1997] states that 
both items have error bars as much as 20%. Furthermore, he states that these random errors will 
skew data to the right when W is in denominator, as seen in Figures 3 and 4. This could tend to 
inflate our oxygen density values when errors are large. 

Error in V|| will increase if vertical winds are present because all measurements assume 
zero vertical winds. A study by Sipler et al. [1995] on nighttime vertical winds above Millstone 
Hill revealed downward movement on the order of 10 m/s with little variability on 
geomagnetically quiet nights (Kp<3) and downward movement on the order of 0 to 50 m/s with 
great variability on geomagnetically active nights (Kp>3), but these values were difficult to 
determine because of error bars ranging from 10 up to 50 m/s. Sipler et al. also noted 
oscillations in the magnitude and direction of the vertical wind velocities. Unfortunately, we 
know little about such oscillations and our data set for this study does not include vertical wind 
information. 

Vertical winds add an additional term to Equation (6) as such: 


= U COS I + sin / — Vfj = W (12) 

where Uz is the vertical wind velocity. Uz will increase our oxygen density calculations because 
it tends to be negative. This error becomes especially significant during geomagnetically 
disturbed periods, when vertical winds may be large and highly variable. We must also consider 
the possible oscillatory nature of the winds when examining the data. 

To determine if any systematic errors, such as listed above, are present in the data, we 
analyze Data Set A and B using the bias formula 



nxO 
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where O is the mean of the MSIS-86 oxygen predictions and n is the number of data points. The 
results are listed Table 2. These indicate a small systematic effect, but the majority of the errors 
are random. 


Conclusion 

This examination of ISR and FPI data provides three significant results: 

1. The data generally supports the MSlS-86 model, although our method displays more 
detail and greater complexity than the MSlS-86 model predicts. Caution should be used, for this 
method also has the possibility for large errors under certain conditions. 

2. A number of nights exhibited trends that were consistently different from MSIS-86 
predictions. Th*ese nights indicate needed refinement in the MSIS-86 model. 

3. Geomagnetic activity does not appear to have a simple one-to-one impact upon the 
average difference between MSIS-86 predictions and the alternative method. Much of the 
impact is disguised either by random errors or the competition between horizontal and vertical 
wind velocities in the end results. 
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In conclusion, employing ISR and FPI data to calculate atomic oxygen density in the 
higher thermosphere provides a ground based method to monitor changes. Unfortunately, the 
method requires simultaneous measurements from large and unique equipment that is available 
at only a few locations in the world. Furthermore, the FPI requires a clear sky. Although limited 
by these factors, this method has the benefit of allowing detailed, long-term study of upper 
thermosphere atomic oxygen density. 
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Table 1. Data Used in this Paper 


Date 

Number of 
Data Points 

Kp .Max 
2 1 00 - 09 00 UT 

Date 

Number of 
Data Points 

Kp Max 
21 00-09 00 UT 

Aug. 19, 1988 

6 

2 

May 27, 1990 

6 

7+ 

Jan. 11, 1989 

7 

2 

Sep. 19, 1990 

12 

4+ 

Feb. 8, 1989 

12 

4 

Sep. 21, 1990 

14 

5+ 

Mar. 1, 1989 

5 

4- 

Sep. 12, 1991 

20 

3 - 

Mar. 9, 1989 

16 

6- 

Oct. 9, 1991 

16 

4+ 

Mar. 20, 1989 

25 

4- 

Mar. 1, 1992 

25 

5 - 

Apr. 13, 1989 

13 

4 

Mar. 2, 1992 

11 

3+ 

Aug. 1, 1989 

4 

2 

May 7, 1992 

28 

3 - 

Aug. 9, 1989 

10 

2- 

Jan. 21, 1993 

38 

3+ 

Oct. 5, 1989 

44 

2+ 

Jan. 26, 1993 

36 

5 

Oct. 6, 1989 

6 

3 - 

Jan. 27, 1993 

35 

3 

Dec. 29, 1989 

52 

5 

Jan. 28, 1993 

11 

2 

Feb. 20, 1990 

52 

5 

Oct. 14, 1993 

3 

3-^- 

Feb. 21, 1990 

79 

5 - 

Mar. 12, 1994 

8 

5+ 

Feb. 22, 1990 

70 

3 

Mar. 21, 1994 

19 

5 

Mar. 19, 1990 

24 

6- 

Apr. 15, 1994 

16 

4 - 

Mar. 22, 1990 

15 

4+ 

Aug. 7, 1994 

5 

1 + 

.Mar. 23, 1990 

5 

4+ 

Aug. 9, 1994 

5 

1 

Mar. 29, 1990 

26 

4- 

Aug. 16, 1994 

4 

3+ 

Apr. 8, 1990 

8 

J - 

Nov. 10, 1994 

J 

5 - 

Apr. 9, 1990 

12 

5-i- 

Dec. 3, 1994 

35 

5 - 

Apr. 10, 1990 

8 

8 - 

Jan. 26, 1995 

26 

1 

May 25, 1990 

22 

4- 

Mar. 3, 1995* 


3+ 


•Data from this day was completely rejected as outliers. 


Table 2. Standard Deviations and Bias 


Set 

Data 

Standard 

Deviation 

S D/Mean 

Bias 

A 

0 

3.1x10''" 

%27.01 



D/W 

7.1x10'“ 

%52.07 

% 18.70 

B 

0 

2.7x10'“ 

%30.24 



D/W 

5.6x10'“ 

%50.69 

%26.1 1 



Data Set A: Aug. 19, 1988 to Mar. 23, 1990 


Figure 1: This scatter plot ot O versus DAV with outliers removed shows the 
wide spread ot data Set A, although a slight linear relationship is apparent. 


Data Set B: Mar. 29, 1990 to Jan. 26, 1995 


Figure 2: This scatter plot of O versus DAV with outliers removed shows the wide 
spread of data Set B. A linear relationship is not as visible as Set A. 








Figure 3; This histogram of DAV with outliers removed shows the 
distribution of Set A. It approximates a log-normal distribution. 



distribution of Set B. It approximates a log-normal distribution. 




25, 


to Mar. 23, 


1990 



Figure 5: This histogram of loa n^u> • u • 

newdistnbutionofSet A. It now approx?ma^e^I 

Fproximates a Gaussian distribution. 
£ ^et B: Mar, 29, 1990 to Jan. 26, 1995 
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36.5 


September 12, 1991 



Figure 7: This chart plots the measured DAV and the predicted O versus the time the data 
was taken. Each O prediction is based upon the estimated altitude the D/W was measured. 
The mea^^urements vary in altitude, as can be seen by the fluctuations in the O predictions, 
but generally range no more than 50 km at about a height of 300 km. 


December 29, 1989 



Universal Time (^OL.rs) 


Figures: Same as Figure 7. Notice the erratic nature of the measured data. 
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Figure 9; Same as Fi 
predicted. 



Figure 10; Same as Fig 
compared to the predicti 
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